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We consider the fluctuations of a time-integrated particle current around an atypical value in a 
generic stochastic Markov process involving classical particles with two-site interaction and hard¬ 
core repulsion on a finite one-dimensional lattice with open boundaries. We address the question 
of which interactions one has to impose on such process to make the atypical value of the current 
typical. It is known that a corresponding effective stochastic Markov process might exist whose 
typical value of the current is equal to the atypical value of the current in the original process within 
a time-translational invariant regime. This effective process has, in principle, non-local transition 
rates. Nevertheless, it turns out that under some conditions the stochastic generator of the effective 
process has the same dynamical rules as the stochastic generator of the original process. We find 
these conditions and show that our approach can be generalized to any time-integrated observable. 
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I. INTRODUCTION 

Rare events and their characterizations are of vital 
importance in different contexts of physics. These phe¬ 
nomena take place on a timescale much larger than the 
timescales characterizing the microscopic dynamics of the 
system. For example nucleation of crystals relies on a 
rare event i.e. the formation of the critical nucleus [l|. 
Protein-folding is also a rare event. In spite of an astro¬ 
nomical number of possible configurations for a protein, 
it folds into a unique native conformation 0. Another 
classical example includes phase transformation for which 
the dynamics might be governed by rare events Q . 

In a general stochastic process the effective interactions 
that induce particular rare events are generally very com¬ 
plicated. In an equilibrium stochastic process the princi¬ 
ple of detailed balance requires that the transition rates 
between a pair of microstates in the canonical ensemble of 
the process satisfy certain relation i.e. the ratio of rates 
for a transition and its time-reverse is given by the Boltz¬ 
mann factor. As a sub-set of this equilibrium ensemble 
one can consider a particular driven ensemble consisting 
of phase-space paths (sometimes called an ensemble of 
trajectories) for which the mean flux of an observable on 
those paths is fixed. The existence of a net flux implies 
that the we are dealing with a constraint driven dynam¬ 
ics. Using the Bayes’ theorem, it has been shown that 
the transition rates of this driven stochastic process with 
a given flux are related to those of the equilibrium sys¬ 
tem [3, @ . Although unphysical transitions in the origi¬ 
nal equilibrium system, which might violate the relevant 
physical laws, will not appear in the driven dynamics; 
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however, the transition rates of the driven system might 
be non-local. The non-equilibrium counterpart to the 
equilibrium detailed balance derived in Si results in a 
set of invariant quantities in the driven system analogues 
to the equilibrium one. This provides us with exact rela¬ 
tions which help us calculate the transition rates in the 
driven system i. The results obtained in S3 can be 
reproduced by maximizing the dynamical entropy in the 
presence of appropriate constraints Q. 

It is long known that in order to study the dynam¬ 
ics of a stochastic process conditioned on atypical val¬ 
ues of a time-integrated observable in the steady state 
of a generic stochastic process system, whether this ob¬ 
servable depends on microstates or transitions between 
a pair of microstates, one can use the concept of bi¬ 
ased ensemble of trajectories @]-[l3- This can be done 
by introducing a biasing field conjugated to the mean 
value of the observable. During a long observation time 
t the ensemble average of a given observable in this bi¬ 
ased ensemble of trajectories might depend on time and 
therefore, the time-translation invariance might be bro¬ 
ken. However, there exists a time interval with 

ti and t 2 being far from the initial time and the final 
time 0 and t respectively, where the time-translation in¬ 
variance is held i.e. the ensemble average of the ob¬ 
servable under investigation in this time interval is in¬ 
dependent of time. It has been shown that being in the 
steady state and during this time-translational invariant 
regime the biased trajectories of the original pro¬ 

cess coincide with unbiased trajectories of an effective (or 
auxiliary) stochastic process [^. Hence the average of 
the observable over the steady-state distribution of the 
effective stochastic process will be equal to its average 
over the biased ensemble of trajectories during the time- 
translational invariant regime. The effective stochastic 
process is a conditioning-free process describing the prob¬ 
lem of conditioning a Markov process on an atypical value 
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of the dynamical observable. The mathematical relation 
between these processes is given by a generalization of 
Doob’s h-transform H- It was shown that the effec¬ 
tive process can be represented as a process satisfying 
various variational principles or a control process opti¬ 
mizing functionals related to the large deviations of the 
conditioning dynamical observable [l5l| . The connection 
between effective interactions and the theory of optimal 
control has also been studied in [I^. The analysis of 
the effective interaction in this way is used, for example, 
in the East model as one of the kinetically constrained 
models consisting of interacting spins in Glass-forming 
systems ll2|. 


The effective stochastic process consists of those in¬ 
teractions one has to impose on the original stochastic 
process to make atypical behavior typical. As in the equi¬ 
librium case explained above, the effective process might 
be unphysical in the sense that its transition rates might 
be non-local Q. This means that the original stochas¬ 
tic process and its corresponding effective stochastic pro¬ 
cess might not share similar features such as the range of 
interactions. The one-dimensional classical Ising chain, 
which exhibits ferromagnetic ordering in its biased en¬ 
semble of trajectories, is an example which reveals this 
feature Q. Similar examples are studied in [13, El- 


A natural question that might arise is that under what 
conditions the corresponding effective stochastic process 
of a stochastic process with conditioned dynamics is 
physical in the sense that, in comparison with the origi¬ 
nal dynamics, no non-local transitions appear in the ef¬ 
fective dynamics. In other words, under what conditions 
imposed on the microscopic reaction rates or for which 
atypical values of the observable, the stochastic genera¬ 
tors of these two processes are exactly the same (up to a 
rescaling of the microscopic dynamical rules). This might 
not be valid for all atypical values of the observable; how¬ 
ever, as we will see one might be able to find at least an 
atypical value of the observable for which the dynamics 
of the original stochastic process and its corresponding 
effective process share the same features at that point. 

In present paper we are going to address the above 
question for a specific class of stochastic Markov pro¬ 
cesses consisting of interacting classical particles on a 
one-dimensional lattice with open boundaries. We as¬ 
sume that the particles are subjected to nearest-neighbor 
interactions in the bulk of the lattice while they can en¬ 
ter or leave the lattice from both the first and the last 
lattice sites. Considering the total reaction-diffusion cur¬ 
rent as a physical observable, we require that the corre¬ 
sponding effective stochastic process consists of exactly 
the same interactions in the bulk and at the boundaries 
of the lattice. In other words, we require that, up to a 
rescaling of the transition rates, the stochastic generator 
of the effective stochastic process is exactly the same as 
the stochastic generator of the original stochastic pro¬ 
cess conditioned on an atypical value of the total particle 
current. We show that, given that there are some con¬ 
straints on the dynamical rules of the original stochastic 


process, there is at most a single atypical value of the 
average current at which this property might be held. 

A couple of examples are given in the present paper. 
In the first example the Asymmetric Simple Exclusion 
Process (ASEP) is considered on an open lattice. In this 
system the particles with hard-core interactions perform 
a continuous-time simple random walk on an open lat¬ 
tice with the possibility of entering or leaving the lattice 
from both the first and the last lattice sites. Consider¬ 
ing a barrier-free hopping of particles between the bulk 
of the lattice and the particle reservoirs with the same 
hopping rates as inside the bulk and assuming that the 
diffusion rates are biased to the right, it turns out that 
the conditions under which the effective and the original 
ASEP share the same features, restrict us to an atypi¬ 
cal value of the particle current which is lower than the 
average particle current in the steady state of the orig¬ 
inal ASEP. On the other hand, it can be seen that the 
effective dynamics is exactly the one for the ASEP but 
with a reversed driving force (i.e. the diffusion rates are 
biased to the left). This phenomenon has already been 
observed in a recent work 0- It has been shown that 
under some constraints the steady state of the effective 
ASEP can be written as a superposition of antishocks. 

In the second example we consider an Asymmet¬ 
ric Kawasaki-Glauber Process (AKGP) on a one¬ 
dimensional lattice with open boundaries 0 . In this 
case the non-zero rates are the death and branching rates 
as well as the hopping rate to the left. It is known that 
stable shocks can develop in the AKGP. Hence, a lin¬ 
ear superposition of them can be used to construct its 
steady-state which consists of a hight-density phase and 
a low-density phase. As we will see by fine tuning the mi¬ 
croscopic reaction rates the stochastic generator of the ef¬ 
fective process can be brought to the form of the stochas¬ 
tic generator of the AKGP. Being in either of the static 
phases, the atypical current at which this phenomenon 
takes place can be lower or higher than the typical value 
of the current in the steady state, depending on the mi¬ 
croscopic reaction rates. 

Finally we will bring the third example in which the 
above mentioned phenomenon can happen for an atyp¬ 
ical value of a non-entropic particle current. While the 
current in the ASEP is entropic and for the AKGP is zero 
(since the steady state is an equilibrium one), interest¬ 
ingly the large deviation function for the current in our 
third example satisfies the Gallavotti-Gohen-like symme¬ 
try [ 2 ^. On the other hand, the atypical value of the 
current at which the above mentioned phenomenon hap¬ 
pens, is exactly equal to the typical value of the current 
in the steady state but with the opposite sign. 

This paper is organized as follows. In section II we 
start with mathematical preliminaries and tools. In sec¬ 
tion III we define the reaction-diffusion current and find 
the conditions under which the stochastic generator of 
the effective process is equivalent with the stochastic gen¬ 
erator of the original process. In section IV we will bring 
three examples to show how our constraints determine 
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the effective dynamics. The generalization is brought in 
section V. The last section is devoted to the outlook and 
conclusion. 

II. MATHEMATICAL TOOLS: A SHORT 
REVIEW 

We start with a stochastic Markov process in 
continuous-time. This is defined through a set of config¬ 
urations denoted by {C} and stochastic transition rates 
ujc^C between these configurations. Considering the 
complete basis vector {IC)}, the probability of finding 
the system in configuration C at time t is given by 
P{C,t) = {C\P{t)) where the ket \P{t)) evolves in time 
according to the following master equation [2l| 

l\p{t)) = n\Pit)) ( 1 ) 

in which the stochastic generator or Hamiltonian P is a, 
square matrix with the following matrix elements 

{C\'H\C') = ujc'^c - ^c,c' ^ wc^c" ■ 

C"#C 

Let us now consider a reaction-diffusion system consisting 
of interacting classical particles on a one-dimensional lat¬ 
tice of length L which is modeled by a stochastic Markov 
process in continuous-time. Being in the steady state, 
we denote the mean (or typical) value of the reaction- 
diffusion current as J*. Let be the number of reac¬ 
tion and diffusion processes which contribute to the total 
reaction-diffusion current of the system up to the time 
t. This quantity is extensive with respect to t and L. 
For a finite L we define the space-time average of the 
total reaction-diffusion current as J = J/{Lt) which is 
a time dependent quantity. During a long-time interval 
t the probability to observe an atypical mean J ^ J* is 
exponentially small in L and t. The large deviation prop¬ 
erty requires P{J) oc exp{—I{J)Lt) where /(J) is called 
the rate function. Now limt_>oo gives the 

cumulant generating function of the current J in which 
= ^je~’^'^PiJ) and that s is called the count¬ 
ing field conjugated to the mean current J Elm. 

We aim to study the dynamics of the above mentioned 
system conditioned on an atypical value of the current 
J. We define Jc^C as an increment for this current 
during transition from configuration C to C . It is known 
that the generating function of J defined above is given 
by = {t\Ps{t)) where (ll| = called the 

summation vector, and that \Ps{t)) should be obtained 
from the following master equation 

= ( 2 ) 

The operator ?i(s) in (I2|) is non-stochastic and called the 
modified Hamiltonian of the system with the following 


matrix elements 

{C\'H{s)\C') = - Sc,C' ^ u)c^c" ■ 

C"/C 

The counting field s can be interpreted as a biasing field 
in the ensemble of dynamical trajectories which is some¬ 
times called the s-ensemble. The role of s in the dynam¬ 
ical ensemble is similar to the parameter /3 (inverse of 
temperature) in the conventional equilibrium canonical 
ensemble. Using this biased ensemble one can study the 
dynamics of system during the observation time t condi¬ 
tioned on a given value of the mean current J. Fixing 
some s ^ 0 correspond to studying those realizations of 
the process in which J fluctuates around some atypical 
mean value This approach is sometimes called the 
grand canonical conditioning which corresponds to con¬ 
structing a canonical ensemble of trajectories [iS[ni[il. 
According to our notation the positive (negative) values 
of the counting field s correspond to the atypical values 
of the current lower (higher) than the typical value of the 
current in the steady state. The time-evolution generator 
or modified Hamiltonian for the conditioned dynamics 
P{s) is a non-stochastic operator which does not con¬ 
serve probability. The sum of unnormalized probabilities 
is called the dynamical partition function of this dynam¬ 
ical ensemble and is given by Z{s,t) = (l|Ps(i))- The 
logarithm of this quantity plays the role of the dynami¬ 
cal free energy of system which determines its dynamical 
phase behavior Q. 

Following the discussion in section HI there is a time- 
translational invariant regime during which one can con¬ 
struct an effective (or auxiliary) stochastic process whose 
unbiased dynamics produces the same value of mean cur¬ 
rent as the conditioned (or biased) ^namics explained 
above does during that time interval l9| . Considering the 
eigenvalue equations for the modified Hamiltonian P{s) 

n{s)\A{s)) = A(s)|A(s)) , 
P(s)(A(s)|=A(s)(A(s)| 

it has been shown that the stochastic generator of this 
effective stochastic process is given by 0 

neffis) = unis)u-^-A*{.s) (3) 

in which 17 is a diagonal matrix with the matrix ele¬ 
ment {C\U\C) = (A*(s)|C') and the asterisk stands for 
the largest eigenvalue and corresponding left and right 
eigenvectors of 77 (s). The off-diagonal matrix elements 
of the operator ?7e//(s) in ([3]) are given by 

(A*(s)|C") 

It is easy to see that for the systems with 
a finite-dimensional configuration space A*(s) = 

limt^oo ln(e“'’‘^)/< [l^. 
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III. EQUIVALENCE OF ORIGINAL AND 
EFFECTIVE DYNAMICS 


in which (j^) is defined as the average local particle cur¬ 
rent from the lattice site fc to fc -I- I and is given by 


In this section we limit ourselves to a family of single¬ 
species reaction-diffusion systems of classical particles 
with nearest-neighbor interactions in the bulk of a one¬ 
dimensional lattice with open boundaries from there the 
particles can enter or leave the lattice. We aim to find 
the conditions under which the effective Hamiltonian 
'Heff{s) of this family is similar to that of the origi¬ 
nal process conditioned on some atypical mean current 
J during its time-translational invariant regime, in the 
sense that the effective Hamiltonian consists of exactly 
the same type of interactions in the bulk and boundaries 
of the lattice. The simplest choice is where 17 in ([3|) is an 
identity matrix. This means that the modified Hamilto¬ 
nian 77 (s) and the effective Hamiltonian 77e//(s) differ 
from each other by a constant which, according to (|5|). is 
the largest eigenvalue of the 77(s). 

For the above mentioned family of stochastic processes 
the Hamiltonian 77 can be written as 

77 = 

+ Efe=i ® h (5) 

-p (g, 7^ 

in which I is a 2 x 2 identity matrix. Introducing the 
basis kets 

"*> = (0' ''‘> = ( 1 ) 

in which 0 and A correspond to a vacancy and an occu¬ 
pied lattice site respectively, the matrix representation of 
h in the basis of {00, %A, A0, AA} and that of £ and iZ 
in the basis of (0, A} are given by 


/ wii LO12 wia UJ14 
^ _ W 21 W 22 W 23 W 24 

W 31 a;32 a;33 a;34 

\ a;4i a;42 a;43 a;44 

The diagonal elements of h are given by oja = 
— ‘j-’ji- As can be seen the parameters a. and 7 (J 

and /3) are the injection and extraction rates of particles 
for the left (right) boundary respectively. 

Let us consider the total reaction-diffusion current as 
the proper dynamical observable. The time-derivative 
of the average local density of particles is related to the 
average particle current through the following continuity 
equation 

= (j’fe-i) - (j’fe) + >5'fe for 1 < k < L ( 6 ) 



(jfc) = [(W21 - W3i)((l - Pfc)(l - Pfc-l-l)) 

— (Wi2 -I-W 42 -I-U;32)((I — Pfc)pfc+l) 

+ (<^ 43 -l-W 13 -|-u;23)(Pfc(l ~ Pfc+i)) (7) 

+ (W24 — W34)(pfePfc+i)](I — (5fe,L)(I — Sk,o) 

+ iPiPk) - S{1 - Pk))Sk,L 

+ (a(l - Pfc+i) - 7(Pfe-si))4,0 

for fc = 0, • • • ,L. Sk is the source term. For the details of 
derivation 0 see Appendix. The average total reaction- 
diffusion current, which includes the contribution of all 
bonds of the lattice, is now given by 

= ( 8 ) 

k=0 


Considering the total reaction-diffusion current defined 
in (|S]) as a dynamical observable, the modified Hamilto¬ 
nian 77 (s) is given by 


77(s) = 

+ (9) 

-p g, 7t(s) 

in which 


/ 

Wii 

wi2e® 

^136 

<j-’l4 \ 


a;2ie“® 

W22 

^236 ^246 ® 


W3ie® 

<^326® 

W33 

W34e® 1 

V 

W41 

a;42e® 

a;43e“® 

^44 / 


£{s) 


—a ye® 
ae“® —7 


,TZ{s) 


—6 /3e ® \ 
(5e® —/3 J 


The increment of the current for each reaction process or 
diffusion process can be understood from ©• Fixing the 
counting field s, corresponding to study of an atypical 
value of the current 77, and trying to find the effective 
Hamiltonian 77e//(s) can be a formidable task. 

Considering {4]) one should note that 77e//(s) can not 
necessarily be written in the two-site interaction form 
though 77 is of the form ([51). Generally speaking, a sim¬ 
ple system might have complex effective interactions. As 
a matter of fact, it has been shown that even for a sys¬ 
tem with short-range interactions the effective interac¬ 
tions might be long-range However, as we 

will see, there might be a value of s = sq at which the 
stochastic Hamiltonian of effective dynamics is similar 
to (O which means it involves nearest-neighbor inter¬ 
actions in the bulk and single-site interactions with the 
reservoirs at the boundaries. 
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We have found that under the following constraints 

^so _ 1^13 + <^23 + ^43 ~ 2c^21 
0^12 + UJ 32 + W 42 — 2 w3i 
a;i3 + a ;23 + ^43 — 0^21 ~ ot 


7 - ^^31 
P - 0 J 21 


UJ12 + a ;32 + 0J42 — W31 — S 
U>24 ~ W21 


( 10 ) 


> 0 


W 34 — a; 3 i 


the effective Hamiltonian TLeffis) has the form of ([51) 
with 


wii 

wi 2 e®“ 

Wise 

Wi4 

a; 2 ie“®“ 

^22 

W23e“®" 

W24e“®“ 

a;3ie®° 

W32e®“ 

^33 

W34e®” 

W 41 

W42e®“ 

W43e“®“ 

W 44 

— 

7 gS 0 




0 

1 

) ’ 


— 

Pe 



- 

-Pe~^° 




heffiso) = 


^e//(so) = 


T^effiso) = 


where the diagonal elements of heff are given by wb = 
— J2j^ii^eff{so))ji which is the stochasticity condition 
for the effective Hamiltonian. The largest eigenvalue of 
H(s) at s = So turns out to be 


A*(so) = (a + {L - l)u} 2 i^ (e - 1) 
+ (, 5 +(L-l)co 3 i)(e««-1) 


( 11 ) 


where its corresponding left eigenvector is given by 
(A(so)l = (1|- One should note that the above left eigen¬ 
vector results in the following exact expression for the 
generating function of the current at s = sq 


^e-soJ) = (l|n„(t)) 


= (l|e‘«^o|P,jo)) 

_ tA*(so) 


( 12 ) 


Depending on the process under investigation, the eigen¬ 
value m might depend linearly on the system size L. 
There are two cases for which the largest eigenvalue 
A*(so) can be independent of the system size. The first 
case is where W 21 = ^31 = 0 while the rest of the reaction 
rules satisfy (fTITl) . In the second case sq = In ^ while 


e “ = 


W 24 _ Wi3 -I- a;23 -b a;43 
W 34 a;i 2 + a;32 -b a;42 
Wi 3 -b UJ 23 + W 43 — a 

7 

P 

W 12 + W 32 + <^42 ~ ^ 


(13) 


IV. EXAMPLES 

In this section a couple of examples are presented to 
show how the conditions obtained in the previous section 
might generate interesting results. In the first example 
we consider the ASEP with open boundaries as explained 
in the introduction. In the bulk of the lattice the particles 
hop to the right and left according to the following rules: 


A% 

0 A 


0 A 
A0 


with the rate 
with the rate 


W 23 = P , 
UJ32 = q 


(14) 


All other reaction rates in the bulk of the lattice are zero. 
The particles are also injected and extracted from the 
boundaries of the lattice with the rates a, 7 , P and 8 as 
explained in the previous section. The constraints (11011 
give 


oSO^P 

q p q 


p 


= 1 . 


(15) 


Let us assume that the density of the particles at the 
left and right boundaries is kept fixed, using two particle 
reservoirs, at the values pi and p 2 respectively. This can 
be done by choosing barrier-free boundary rates defined 
as [ 1 ^ 

a=ppi, 7 = g(l-pi), P=p{l-p 2 ), S = qp 2 . 

In this case the only constraint which remains will be 


q 

and the eigenvalue is given by 

A*(so) = -{p-q){pi - P2) ■ 


(16) 


(17) 


By substituting sq in ^(s) it is easy to see that the effec¬ 
tive Hamiltonian can be obtained from the Hamiltonian 
of the original ASEP by exchanging p and q (or reversal 
p q oi the particle hopping rates). This observation 
has an interesting consequence. Given that one chooses 
barrier-free boundary rates, it is known that the steady 
state of the ASEP can be written as a linear superposi¬ 
tion of Bernoulli measures with a step-function structure 
provided that the following constraint is satisfied [l^ 

P2(l - Pi) ^ P 
Pi(l - P 2 ) q ' 

Hence, if one chooses p > q then the constraint requires 
Pi < P 2 and this is what we call a shock structure. Now, 
since the effective Hamiltonian is exactly the same as 
the original Hamiltonian but with reversed hopping rates 
one can conclude that the steady state of the effective 
Hamiltonian can also be written as a superposition of 
Bernoulli measures with a step-function structure. In 
this case one should have 

P2(l - Pi) ^ q 

Pi{l-P2) P 


> 0 . 
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and since p > q then pi > p 2 - In comparison to the 
definition of a shock structure this is called an antishock. 
One should note that since ior p > q we have sq > 0, 
then atypical value of the current is always lower than 
the typical value in the steady state. This has already 
been observed and discussed with more detail in 0- 
In the second example we consider an asymmetric 
Kawasaki-Glauber process which contains the following 
reaction rules in the bulk of the lattice: 


%A 00 
0 A ^ A% 
%A AA 
A0 00 
A0 AA 


with the rate 
with the rate 
with the rate 
with the rate 
with the rate 


Wl 2 , 

^32 , 

UJ42 , (18) 

Wl3 , 

W 43 • 


The only non-zero boundary rates a and /3 define the in¬ 
jection and extraction of the particles at the left and right 
boundaries of the lattice respectively. The constraints 
for this process are 



-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 

S 


FIG. 1 : The plot of the numerically calculated largest eigen¬ 
value of the modified Hamiltonian. The dotted curve is 
Aturi^) for i^>i2 = l,a;i3 = '2,uj21 = 0.8 for a lattice of 
length L = 6. The solid line and the dashed line correspond 
to A|„t(s) and respectively. The vertical line is 

So = ln(a;i3/tJi2) = 0 . 69 . See inside the text for more infor¬ 
mation. 


e 


So 


Wi 3 + UJ43 a 

- ;-;- , a = ji = UJl3 + LO43 ■ 

W 12 + ^42 -|- a;32 


It is known that the steady state of this process, without 
any constraints on the microscopic reaction rates, can be 
written as a superposition of stable Bernoulli shock mea¬ 
sures [ 1 ^. It has also been shown that the microscopic 
position of each shock performs a biased random walk 
on the lattice. Now, following our discussion in the first 
example, we conclude that the steady state of the effec¬ 
tive dynamics can be written in terms of superposition 
of Bernoulli shock measures (and not antishocks). Note 
that in the steady state of the original process the system 
undergoes a static phase transition between a low-density 
and a high-density phase depending on the values of a;i 3 
and a; 43 . It is worth mentioning that, being in either of 
these static phases, the dynamics can be either condi¬ 
tioned on a lower than typical or a higher than typical 
value of the total average current. 

The time-integrated currents are generally either en- 
tropic which satisfy the Gallavotti-Cohen symmetry such 
as the one studied in or non-entropic which sat¬ 
isfy the Gallavotti-Cohen-like symmetry such as the one 
studied in [^ . In the third example we introduce 
a non-entropic reaction-diffusion current which satisfies 
the Gallavotti-Gohen-like symmetry with the mentioned 
property. It turns out that under some constraints the 
value of the conjugated field sq can be located on the 
symmetry point corresponding to s = 0 which means we 
have A*(0) = A*(sq 7^ 0) = 0. This indicates that the 
absolute values of the atypical current and the typical 
current are equal; however, they flow in opposite direc¬ 
tions. 

Our third example consists of birth and death pro¬ 
cesses in the bulk of the lattice with the rates (a; 2 i,a; 3 i) 
and (wi2, W13) respectively which can be demonstrated as 


follows: 


00 - 

0 A 

with the rate 

W21 

00 - 

A A0 

with the rate 

W31 

%A - 

-A 00 

with the rate 

^12 

A0 - 

-A 00 

with the rate 

^13 


(19) 


The particles are allowed to enter or leave the lattice from 
both boundaries; however, the boundary rates a and 5 
are assumed to satisfy the following constraints 


a = UJ 21 , 6 = a; 3 i . 
Now the constraints (fTOl) lead us to 

a; 2 i a;i3 

So = In-= In-, 

wsi a;i 2 

7 = ^12 — W31, 

l3 = 0Ji3- a;2i ■ 


The largest eigenvalue of the modified Hamiltonian for 
the entropy production A*„((s) and that of the total 
diffusion-reaction current are numerically cal¬ 

culated and plotted in FIGHJ We have also plotted 
A*„^(£'s) where E = sq. As can be seen does not 

lie on A*^^{Es) and therefore the total diffusion-reaction 
current is non-entropic [^ . 


V. GENERALIZATIONS 

The above discussion can be generalized to any arbi¬ 
trary time-integrated observable (which is not necessar¬ 
ily the particle current) in a continuous-time stochastic 
Markov process with a stochastic generator of type ([5]) 
and a finite configuration space. These observables can 
be fluxes or currents which depend on the transitions 
between configurations or microstates, such as the one 
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we explained in this paper. Alternatively we can con¬ 
sider those time-integrated observables that might have 
merely a spatial nature such as dynamical activity Q or 
energy 0. 

For the fluxes or currents which are defined on the ba¬ 
sis of transitions between configurations, we consider the 
increment Oc^c whenever the system jumps from C to 
C along a spatio-temporal trajectory. For the dynamical 
activity one has Oc^c = 1 for all C and C {C ^ C) 
while for the entropy production the increment will be 
Oc^c = dll- We have already de¬ 

fined the increments for a global reaction-diffusion cur¬ 
rent in section lllll These increments, as we saw, affect the 
non-diagonal elements of the modified Hamiltonian. In 
contrast, for those time-dependent observables which are 
defined along a spatio-temporal trajectory and depend 
on the visited microstates, only the diagonal elements of 
the modified Hamiltonian, depending on the observable, 
are changed. 

In either of these two cases we start with constructing 
the modified Hamiltonian 77(s) for the observable under 
investigation. Let us denote the sum of the matrix ele¬ 
ments of ith column of h(s) as hi for i = 1, • • • ,4. For 
7^(s) and C{s) they will be denoted by and h respec¬ 
tively for 1 = 1,2. It can be shown that the summation 
vector (1| is the left eigenvector of the modified Hamil¬ 
tonian given that 

hs — hi = hi — h2 = h — h = r2 — ri, h/i = hi . ( 20 ) 

At the same time the eigenvalue of the modified Hamil¬ 
tonian associated with that left eigenvector is given by 

A*(s) = li ri -\- {^L — l)/ii (21) 

in which L is the size of the lattice. The equations (1201) 
determine the value(s) of the conjugated field sq and 
also the probable constraints on the microscopic reaction 
rates under which the original and the effective dynam¬ 
ics are equivalent in the sense of what was explained in 
section Ell 


VI. CONCLUDING REMARKS 

In order to investigate the dynamics of a generic 
stochastic Markov process conditioned on an atypi¬ 
cal value of an integrated current during its time- 
translational invariance regime, one can modify its 
stochastic generator to build an effective (or auxiliary) 
stochastic generator for which the typical value of the in¬ 
tegrated current in the steady state is equal to the atypi¬ 
cal value of the integrated current in the original process. 
However, one realizes that the resulting effective process 
might be unphysical in the sense that it might contain 
non-local transitions. In this paper we have shown that 
under some constraints on the microscopic reaction rates, 
the stochastic generator of the effective stochastic pro¬ 
cess can posses exactly the same dynamical rules as the 


original process does, at least for a specific value of the 
current under investigation. We have also shown that, 
depending on the process, this current might be entropic 
or non-entropic. Possible generalizations have also been 
discussed. Our approach might not be the only pos¬ 
sible way to construct such effective stochastic process 
who shares identical features with the original stochastic 
process. It would be of great interest if one could find 
the general conditions under which the effective process 
would be physical in the sense that it only contains local 
transitions. On the other hand, we only considered the 
reaction-diffusion processes with nearest-neighbor inter¬ 
actions on open lattices. It would be interesting to in¬ 
vestigate the processes with long-range interactions, not 
only on an open lattice but under periodic boundary con¬ 
ditions. 


Appendix: Derivation of the particle current 
formula (Q 


The time evolution of the average local particle density 
{pk){t) at the lattice site k at time t is given by 




.fc + Jk->^k+l + Jk-l,k ~ Jk,k+l (^-1) 


for fc = 1 , 2 ,..., L where J^_i^i^ and are the av¬ 

erage input current into the lattice site fc, in the result of 
reaction with the lattice sites k — 1 and k + 1 respectively. 
'^kk+i '^k-i k average diffusion cur¬ 

rent from the lattice site fc to fc -|- 1 and from fc — I to 
k respectively. Note that = 0 while 

J(Pi = i^^i 0 which give the particle exchange with 
particle reservoirs at the boundaries. These quantities 
are given by 


jR 

a; 


k — l—¥k 


J, 


R 


k^ —fc-j-l 


= (W21 + W4l)((l — Pfe_l)(I — Pfe)) 

— Wl2((I - Pk-l)pk) + W43(Pfc-l(l - Pk)) 

— (Wi4 -I-W34)(pfc_ipfc) (I — (jfc,l) (A.2) 

= (W3i -k W4l)((l - Pfc)(l - Pfc+l)) 

+ W42((I — Pk)Pk+l) — Wl3(Pfc(l “ Pk+l)) 

— (Wi4 + W24)(PfcPfc+l) (1 — <5fc,L) (A.3) 


for fe = I, 2,..., L and 


Jk,k +1 = W23(Pfe(l - Pfc+i))(I - 4 .l)( 1 - 4,o) 

— W32((l — pfe)pfe+i)(l — (5fe,L)(l — 5kfi) 

+ (/?(Pfe) - - Pfc))4.L 

-I- (a(I - Pfc+i) - 7(pfc-i-i))4,o (A.4) 

for k = 0, • • • ,L. The average local density of particles 
is related to the average particle current through the fol¬ 
lowing continuity equation Comparing (lA.ll) and dH]) 








one finds the following relations for the average particle 
current (j^) and the source term Sk 


{jk) — Jk^k+l + '^k,k+l J 


Sk = + Jf" 


k—¥k-\-l 
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